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We use first principles density functional theory to investigate the softening of polar phonon 
modes in rutile TiC^ under tensile (llO)-oriented strain. We show that the system becomes unstable 
against a ferroelectric distortion with polarization along (110) for experimentally accessible strain 
values. The resulting polarization, estimated from the Born effective charges, even exceeds the 
bulk polarization of BaTiOs- Our calculations demonstrate the different strain dependence of polar 
modes polarized along (110) and (001) directions, and we discuss the possibility of strain engineering 
the polarization direction, and the resulting dielectric and piezoelectric response, in thin films of 
Ti02 grown on suitable substrates. 
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Rutile Ti02 is an incipient ferroelectric materialf^ 
and it has been shown via ab initio simulations that a 
ferroelectric transition can be induced by either negative 
pressure or (OOl)-oriented strainP^ An experimental re- 
alization of ferroelectric Ti02 would have great techno- 
logical impact, since this material is cheap and optimized 
processing techniques are readily available. Ti02 is cur- 
rently used for a large variety of applications such as solar 
cells or pigments. Besides, due to its high dielectric con- 
stant (e = 251 at low temperature^), Ti02 is also used in 
dielectric devices. Hence, the possibility to tune the elec- 
tric permittivity with or without inducing a ferroelectric 
transition would be of great technological importance. 

Most applications are based on the (110) surface, which 
can be strained in either (110) or (001) direction by 
growth on a suitable substrate. Until now, only the in- 
fluence of isotropic and (OOl)-oriented strain on a fer- 
roelectric state with polar ization along (001) has been 
investigated for bulk rutile pSMU The resulting ferroelec- 
tricity can be attributed to the only polar phonon mode 
polarized in (001) direction, A 2n , which softens if the Ti- 
O distances are enlarged. Interestingly, it was predicted 
recently that the A 2u mode is insensitive to (110) strain 
whereas an acoustic phonon mode softens under these 
strain conditions P Apart from A 2n , rutile possesses ad- 
ditional low energy phonon modes which are polarized 
along (110), and for which a detailed understanding is 
still missing. 

Here, we investigate the role of (llO)-oriented strain 
on the ferroelectric properties of rutile Ti02- The main 
goal of our work is a systematic investigation of the in- 
fluence of (110) strain on polar modes in (001) and (110) 
direction. We show that the paraelectric state of rutile is 
destabilized under tensile strain and a polarization along 
(110) is induced. Furthermore, a polarization in (001) di- 
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rection can be obtained if the ionic shifts along the (110) 
direction are suppressed. 

All results presented here are obtained from self- 
consistent density functional theory calculations employ- 
ing the "Vienna Ab Initio Simulation Package" (VASP) 9 
with a plane wave basis, projector augmented wave 
potentials,^ and the generalized gradient approxima- 
tion (GGA) in the formulation of Perdew, Burke and 
Ernzerhof. 11 For Ti (O) atoms the 4s3d (2s2p) electrons 
are treated as valence. Ionic positions are optimized until 
all forces are smaller than 0.01 eV/A and a \[2 x \[2 x 1 
supercell containing 4 Ti02 units, see Fig. [I] is used. 
Born effective charges have been calculated using den- 
sity functional perturbation theory.^ 

The rutile crystal structure of Ti02 is tetragonal with 
space group symmetry P^/mnm (see the primitive unit 
cell inserted in Fig. [l^). Each Ti-atom is 6- fold coor- 
dinated by O-atoms, with 4 short equatorial Ti-0 bonds 
and two longer apical bonds, see Fig.[l]o. The exact Ti-0 
distances are determined by the lattice constant a, the 
tetragonal c/a ratio, and the internal structural param- 
eter u. From our relaxations we obtain a = 4.6640 A 
(4.587 A), c/a = 0.6366 (0.644), and u = 0.3047 (0.305). 
For comparison, experimental results from Ref. 13 , using 
neutron diffraction at 15 K, are given in brackets. 

As mentioned in the introduction, one polar optical 
A 2u mode with polarization along (001) exists in rutile. 
The corresponding ferroelectric state can be stabilized by 
increasing the equatorial Ti-0 bond length, since this re- 
duces the short range repulsion between the two ionsP^ 
Furthermore, there are three pairs of doubly-degenerate 
polar optical phonon modes of E u symmetry, which cor- 
respond to atomic displacements within the (HO)-(llO) 
plane. 1 In Ref. 3 it has been shown that a uniform lattice 
expansion causes a softening of the energetically lowest 
E u mode, although this effect is smaller than the corre- 
sponding change for the A2 n mode. 

The application of (110) strain reduces the symmetry 
to the orthorhombic space group Cmmm, with two in- 
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FIG. 1. Atomic structure of bulk rutile in the \f2 x \f2 x 1 supercell. Solid arrows indicate the atomic shifts for the three (110) 
polarized i^-modes used in our analysis, sorted in ascending order of their energies, (a) Bzu,\ mode, indices mark inequivalent 
Ti/O sublattices under (110) strain, and the dotted circle marks the O atom which has been shifted previous to the structural 
optimization; (b) B^u^i mode; (c) B^ u ^ mode. Dotted lines and arrows mark the primitive unit cell, the lattice constants a 
and c, and the internal structural parameter u. 



equivalent Ti (Tii,Ti2) and O (Oi, O2) sublattices. The 
imposed strain changes the apical Tii-02 and the equa- 
torial Ti 2 -02 bond distance, whereas the equatorial Tii- 
Oi and the apical Ti2-Oi bonds are not modified, see 
Fig. [TJa). As a result, the twofold degeneracies among 
the six optical E u modes are lifted and they split into 
three £?2u and three B^ u modes with atomic shifts along 
(110) and (110), respectively. 

For further analysis of the atomic displacements it is 
convenient to decompose the 3 x N dimensional displace- 
ment vector Y into a linear combination of three mode 
patterns ~Bs u ,i and an additional relaxation R, 

Y = A • B 3 „,i + B • B 3n , 2 + C • B 3n , 3 + R. (1) 

Here, N is the number of atoms within the supercell. We 
use a set of B 3ui modes with displacement patterns as 
sketched in Fig. [I] and the same magnitude of displace- 
ments on all sublattices. These three modes, together 
with the corresponding acoustic mode, form a complete 
basis for representing the B% u eigenmodes of the strained 
rutile structure. Since we are only interested in the three 
optical modes, we do not consider the acoustic mode in 
our analysis. We note that even though the displacement 
amplitudes of the various sublattices in the exact phonon 
eigenvectors cannot be fully determined by symmetry ar- 
guments alone, the three mode vectors depicted on Fig.[l] 
are indeed relatively close to the actual E u eigenvectors 
of unstrained Ti02- 

To explore the energy landscape of (llO)-strained TiC>2 
as function of the B% u optical modes, we perform to- 
tal energy calculations with fixed atomic positions, and 
include different amounts of the three distortion modes 
Bsu,i- The system is softest along B^ u ,i^ but a combina- 
tion of Bsu,i and Bs u ,2 is required to produce a stable 
ferroelectric energy minimum, see Fig. [2] This energy 
minimum is shifted to larger amplitudes for increasing 
strain (compare Figs. and [2)3), but no major differ- 
ences appear for higher strain values. In the following we 



focus our discussion on the case of 2 % strain. 

In order to locate an approximate energy minimum, we 
first identify the most favorable B 3u ,2 amplitude as func- 
tion of £?3u,i, and then vary the amplitude of Bs Uy s for 
the resulting optimal linear combinations of B^ u ,i and 
£>3u,2- The so-obtained energy landscape is shown in 
Fig.[2]3. From this we derive an approximate energy min- 
imum at Yo corresponding to A = 0.096 A, B = 0.035 A, 
C = 0.017 A, andR = 0. 

The results of our full structural relaxations for different 
magnitudes of (110) strain are listed in Table [l[ The sym- 
metry of the system has been reduced to Pmm2 prior to 
the relaxation by shifting one of the 2 atoms 0.002 A 
along (110), see Fig. [IJa). 

The atomic relaxations confirm that for tensile (110) 
strain the paraelectric state is destabilized by a polar 
shift of the sublattices against each other in (110) di- 
rection. This shift, which is well approximated by Yq 
(see Table [i]), leads to alternating expansion/contraction 
of the short Ti2-02 bonds. However, since a further in- 
crease of this bond length (besides the increase resulting 
from the overall (110) strain) is unfavorable, additional 
atomic displacements corresponding to an M point zone 
boundary mode lead to different (110) shifts of the two 
O2 atoms. Furthermore, a small relaxation of the Oi 
atoms in (110) direction (not included in Table [i]) is also 
superimposed to Yo- As a result, the Tii-Oi bond length 
is approximately conserved, even though these ions are 
shifted against each other in (110) direction. This lowers 
the energy for the (110) shift and thus the amplitudes of 
the (110) shift for these sublattice are slightly increased 
compared to Yo- 

The resulting spontaneous polarization can be esti- 
mated by 

P (iio) = y^2x(i)-Z^ 110)ti . (2) 
Here, the summation is over all atoms in the unit cell with 
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FIG. 2. Energy landscapes for static displacements along the 
Bzu,% mode vectors corresponding to 2 % strain (a),(c) and 
5% strain (b). Amplitudes of the modes are given in A and 
energies in meV/Ti02 unit relative to the unrelaxed strained 
configuration. 



volume V, x(i) represents the shift of atom i in (110) di- 
rection, and ^* 110 ) i the corresponding principal value of 
the Born effective charge tensor for undistorted bulk. As 
the Born charges in Ti02 are quite stable with respect 
to small modifications of the Ti-0 distances,^ this ap- 
proximation is sufficient for a qualitative discussion. We 
obtain Z* no) = 7.59|e| (— 5.14|e|) for Ti x (0 2 ), which has 
the apical bond aligned along (110), whereas we obtain 
Z (*no) = 5 - 35 l e l (- L32 l e l) for Ti 2 (Oi) which is in quali- 
tative agreement to Z* 110 ^ = 7.34|e| (— 4.98|e|) along the 
aptical bond, respectively Z* 110 ^ = 5.34 1 e | (— 1.36|e|) for 
the perpendicular direction obtained by Lee et al. ffl 



TABLE I. Amplitudes x (in A) of atomic shifts in (110) di- 
rection under tensile strain. The energy differences AE rel- 
ative to the paraelectric state at the same strain are given 
in meV/TiCb unit, and the resulting polarization P s in 
/iC/cm -1 . 



Strain 


*(Tii) 


x(Ti 2 ) 


x(Oi) 


*(0 2 ) 


AE P s 


1.0^ 


0.09 


0.04 


-0.02 


-0.04 


4.8 33 


1.02 


0.10 


0.04 


-0.03 


-0.03/-0.05 


6.0 35 


1.03 


0.13 


0.05 


-0.03 


-0.04/-0.08 


15.2 46 


1.04 


0.16 


0.05 


-0.03 


-0.04/-0.10 


29.8 54 


1.05 


0.19 


0.05 


-0.03 


-0.05/-0.13 49.3 62 



a Shift along Yq without further atomic relaxation. 



With Eq. [2] we calculate a sizeable polarization, which 
is already in the technical relevant range for a strain of 
2 %, see Table |lj The resulting ferroelectric state is 6 
meV/formula unit lower in energy than the paraelectric 
state. If one assumes that the ferroelectricity is related 
to exactly one degree of freedom at the T point this cor- 
respond to a thermal energy of about 139 K. However, 
since the depth of the energy well increases strongly with 
the imposed strain, a ferroelectric phase should be exper- 
imentally accessible for about 5 % tensile strain as our 
rough estimation already results in a thermal energy of 
1146 K in this case. 

For compressive (110) strain, the reduction of the Ti- 
O distances induces a large short range repulsion, and 
a ferroelectric distortion is unfavorable. Instead, O2/O1 
atoms increase the (110)/(-110) component of their equa- 
torial bond, due to the larger compressibility of the apical 
bond. This is in agreement with the hardening of the fer- 
roelectric A2n mode under compression.^ 

We note that we could not reproduce the destabiliza- 
tion of the system along an acoustic phonon mode un- 
der (110) strain found in Ref. 8, for which an energy 
gain of several meV has been predicted, even if we use a 
2 • y/2 x 2 • y/2 x 4 supercell which is commensurable with 
the displacement pattern of this mode. Possible expla- 
nations for this discrepancy are the different potentials 
used and the resulting difference in lattice parameters. 

Finally, we discuss the relationship between the polar 
modes with polarization along (110) and (001), respec- 
tively. Since the equatorial Ti-0 bond length increases 
under tensile (110) strain, the short range repulsion be- 
tween these ions is reduced, and a polar shift along (001) 
should become likely. Nevertheless, no softening of the 
A2 U -mode under (110) strain has been observed in pre- 
vious ab initio investigations P In contrast, our calcula- 
tions show that the polar mode in (001) direction can be 
stabilized if the displacements along the B% u modes are 
suppressed, see Fig. U 14 If the ions are not allowed to 
displace along (110), a characteristic double well poten- 
tial appears for tensile (110) strain, and the well depth 
increases with increasing strain. However, a coupling 
between soft modes corresponding to (110) and (001) 
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FIG. 3. (a) Total energy change for static atomic displacements along the A2n-mode under uniaxial (110) strain for atomic 
positions relaxed with imposed Cmmm/ Pmm2 symmetry. Here, only the second symmetry class allows for polarization in 
(110) direction, (b) Energy landscape for mixed displacements along Yo and Aiu for 2 % tensile strain in (110) direction. 



polarization exists, similar to perovskites such as e.g. 
SrTiOsP^ As a result, the ferroelectric mode polarized 
in (110) direction disables polar distortions in (001) di- 
rection. If the strain increases from 2% to 5%, and the 
relative ionic shifts in (110) direction increase, more en- 
ergy is needed to shift the sublattices relative to each 
other along (001). Fig.[3]3 shows the energy landscape as 
function of polar shifts along (001), A 2n , and along the 
optimized mode vector Yo in (110) direction. 

While for pure (llO)-oriented strain the energy gain 
corresponding to the ^n-mode is one order of magnitude 
smaller than the energy gain corresponding to B^ U1 a fur- 
ther softening occurs under tensile (001) strain,^ which 
is stronger for A^ u than for B^ u . The exact location 
and depth of the energy minimum within the U A2 U -Bs u 
plane" can therefore be adjusted by straining in both 
(110) and (001) directions, which can be achieved in 
thin films by choosing an appropriate substrate. This 
opens up the exciting possibility of "strain engineering" 
the polarization direction in Ti02, and the resulting di- 
electric and piezoelectri c resp onse, similar to the case of 
perovskite ferroelectrics P^USl 

In summary we have shown that a ferroelectric state 
can be stabilized in rutile under tensile (llO)-oriented 
strain. The main contribution of the displacement pat- 
tern can be attributed to the E u phonon modes of undis- 
torted rutile. For 5% strain a large polarization of 



63/iC/cm -1 emerges in (110) direction. The depth of 
the corresponding energy well is large in comparison to 
thermal fluctuations at low temperatures, which suggests 
that the ferroelectric state can indeed be observed ex- 
perimentally. We note that while the well-known small 
overestimation of unit cell volumes within GGA may also 
slightly overestimate the tendency towards ferroelectric 
distortions, qualitative trends are nevertheless described 
properly and agree well with other calculations! 4 * 6 * 17 -^ 

Furthermore, we have demonstrated that a ferroelec- 
tric state in (001) direction can also be stabilized under 
tensile (110) strain, if the polar shifts along the (110) di- 
rection are suppressed. This indicates a coupling between 
optical modes polarized in (001) and (110) directions, 
which allows to shift the relative energies of ferroelectric 
states with polarization in (110) and (001) directions un- 
der different strain conditions. Future investigation are 
necessary in order to investigate the influence of addi- 
tional (001) strain on the ferroelectric modes in more 
detail. In addition, the effect of surface-induced atomic 
relaxations at the real (110) surface on the polar phonon 
modes is important and will be addressed in an upcoming 
investigation.^ 
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